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Abstract. We review the conformal field theory approach to entanglement 
entropy in 1+1 dimensions. We show how to apply these methods to the 
calculation of the entanglement entropy of a single interval, and the generalization 
to different situations such as finite size, systems with boundaries, and the case 
of several disjoint intervals. We discuss the behaviour away from the critical 
point and the spectrum of the reduced density matrix. Quantum quenches, as 
paradigms of non-equilibrium situations, are also considered. 
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1. Introduction 

Entanglement is one of the most fundamental and fascinating features of quantum 
mechanics, yet in some ways the most mysterious: performing a local measure may 
instantaneously affect the outcome of local measurements far away. This phenomenon 
has been the basis for the development of such new branches of research as quantum 
information and communication. A very recent and rich field of research concerns the 
understanding of the role of entanglement in many-body systems. 

Beside its own fundamental theoretical interest, a principal reason for the success 
of the entanglement entropy as an entanglement measure in extended quantum systems 
is surely its universal scaling at one-dimensional (ID) conformal critical points. The 
equation [1, 2, 3] 

S A = C -ln- + c' 1 , (1) 

O (X 
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has become one of the most ubiquitous formulas in the last five years' literature, 
appearing in fields as apparently unrelated as quantum information, condensed matter, 
and high energy physics. The reasons for this prominence are clear: it is a single 
quantity, easily measurable in numerical simulations, that at the same time gives the 
location of the critical point and one of its most important universal signatures, the 
central charge c of the underlying conformal field theory (CFT). 

The aim of this review is to give a self-contained presentation of most of the 
results for the entanglement entropy that can be obtained by means of CFT. Some 
other important features of the entanglement in extended quantum systems will only 
be considered on passing. For a comprehensive treatment of these aspects we refer the 
reader to the other reviews in this special issue [4] and to the already existing ones 
[5, 6]. 

The plan of this article is as follows. In Sec. 2 we review the CFT (and more 
generally quantum field theory) approach to entanglement entropy based on the 
replica trick, the mapping to the partition function on Ricmann surfaces, and the 
introduction of twist fields. In Sec. 3 we apply these methods to the calculation of the 
entanglement entropy of a single interval, showing in particular how to obtain Eq. (1), 
and generalizing it to several different situations, like finite size, finite temperature, 
systems with boundaries and defects. In Sec. 4 we consider the case of several disjoint 
intervals, with particular attention to the case of two intervals, where several results are 
now available. In Sec. 5 we consider massive perturbations to the conformal behaviour 
in the regime when the mass is small and the systems still retain signatures of the 
close conformal critical point. In Sec. 6 we derive the consequences of the conformal 
scaling for the full spectrum of the reduced density matrix. In Sec. 7 we discuss the 
CFT approach to non-equilibrium situations known as quantum quenches. Finally in 
Sec. 8 we report on an interesting proposal of Klich and Levitov [7] to measure the 
entanglement entropy in real experiments. 



2. Entanglement, replicas, Riemann surfaces, twist fields and all that 

2.1. Basic definitions 

Let p be the density matrix of a system, which we take to be in the pure quantum 
state \^f), so that p — \^)(^f\. Let the Hilbcrt space be written as a direct product 
TL = Ha <8> H.B- A's reduced density matrix is pa = Tr Bp- The entanglement entropy 
is the corresponding von Neumann entropy 

S A = ~Tr p A lnp A , (2) 

and analogously for Sb- When p corresponds to a pure quantum state Sa = Sb- For 
future use, we also define the Renyi entropies 

S^ = T L-]nTrp n A , (3) 

that also are characterized by S A n ^ = S^ whenever p corresponds to a pure quantum 

state. From these definitions S A = lim n _>i S A n \ 

When a system is in a mixed state the entanglement entropy is not any longer a 
good measure of entanglement since it clearly mixes quantum and classical correlations 
(e.g. in an high temperature mixed state, it must reproduce the extensive result for 
the thermal entropy that has nothing to do with entanglement.) This is also evident 
from the fact that Sa is no longer equal to Sb- A quantity that is easily constructed 
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from the knowledge of Sa and Sb is the so called mutual information, defined from 
the Renyi entropy as 



be an entanglement measure (see the review by Amico and Fazio in this volume [8] for 
a discussion of various entanglement measures), but it has the fundamental property 
of satisfying the area law [9, 6] even at finite temperature [10]. 

2.2. A replica approach 

When dealing with a statistical model with a finite number of degrees of freedom, the 
most direct way to obtain the entanglement entropy is to construct the reduced density 
matrix (or at least its eigenvalues Xi as in DMRG) exactly or numerically and then, 
by brute force or by analytic methods, calculate the sum Sa — — X) ^» m Aj. Several 
examples of how this can be worked out, even analytically, for the simplest models are 
reported in other reviews in this special issue. However, calculating the full reduced 
density matrix for a generic interacting quantum field theory remains a daunting 
challenge from an operational point of view, and so here we will take a different route, 
following our previous papers [3, 11, 12, 13]. The approach is reminiscent of the 
"replica trick" in disordered systems and was present in an embryonic form also in the 
early paper by Holzhey, Larsen, and Wilczek [1] . Let us start by considering a lattice 
model. The eigenvalues of the reduced density matrix Xi lie in the interval [0, 1] and 
^2 Xi = 1. Thus, for any n > 1 (even not integer), the sum Tr p\ = ^\ A™ is absolutely 
convergent and therefore analytic for all Ren > 1. The derivative wrt n therefore also 
exists and is analytic in the region. Moreover, if the entropy Sa — — X)i 1°S ^ s 
finite, the limit as n —* 1 + of the first derivative converges to this value. Thus if we 
are able to calculate Tr p n A for any n > 1 we have the entanglement entropy as 



However, calculating Tr p n A for a generic real n in a quantum field theory is still 
an hopeless task. And here the "replica trick" enters: we compute Tr p A only for 
positive integral n and then analytically continue it to a general complex value. We 
will see that the calculation of Tr p A reduces to that of a partition function on a 
complicated Riemann surface (or equivalently to the correlation function of specific 
twist fields) that is analytically achievable in a quantum field theory. The problem is 
then moved to the existence and the uniqueness of a proper analytic continuation to 
extract the entanglement entropy (as in any approach based on replicas). In some cases 
this is straightforward, in others difficult, and in several others beyond our present 
understanding. 

2. 3. Path integral formulation and Riemann surfaces 

To begin with a well-defined problem, let us consider a lattice quantum theory in 
one space and one time dimension (the generalization to higher spatial dimensions is 
straightforward). The lattice sites are labelled by a discrete variable x and the lattice 
spacing is a. The domain of x can be finite, i.e. some interval of length L, semi-infinite, 
or infinite. Time is continuous. A complete set of local commuting observables will 
be denoted by {(f> x }, and their eigenvalues and corresponding eigenstates by {<ft x } and 



r(«) _ c(«) i c(«) c( : 

1 A:B ^ A ' ,:, A 




that is by definition symmetric in 




(5) 
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Figure 1. From density matrix to reduced density matrix. Left: Path integral 
representation of p(<j>\<j>'). Center: The partition function Z is obtained by sewing 
together the edges along r = and r = f) to form a cylinder of circumference (3. 
Right: The reduced density matrix pj\ is obtained by sewing together only those 
points x which are not in A. 



| {4> x }) respectively. (For a bosonic lattice field theory, these will be the fundamental 
bosonic fields of the theory; for a spin model some particular component of the local 
spin.) The states ^K^z}) = | risi^x}) form a basis. The dynamics of the theory is 
described by the Hamiltonian H . The elements of the density matrix p in a thermal 
state at inverse temperature (3 are 

p({<p x }\{<p' x ,}) = (Q{<MIpI n^'» = z{(3y 1 (Y{{^\^ H \ ' ( 6 ) 

x x' x x' 

where Z((3) = Tre~ 0H is the partition function. This may be written as a path 
integral on the imaginary time interval (0,(3): 

p{{4>x}MA) = Z- 1 J [&f>(y, t)\ 11 5{4>{y, 0) -#,,)![ Wiv, 0) ~ M ^ Se , (7) 

X 1 X 

where the euclidean action is Se = f Ldr, with L the euclidean lagrangian. This 
is illustrated in the left panel of Fig. 1. Here the rows and columns of the reduced 
density matrix are labelled by the values of the fields at r = 0, (3. 

The normalization factor Z is the partition function, and ensures that Trp = 1. 
It is found by setting {4> x } = {4>'x\ an d integrating over these variables. In the path 
integral, this has the effect of sewing together the edges along r = and r = (3 to 
form a cylinder of circumference (3 as depicted in the center of Fig. 1. 

Now let A be a subsystem consisting of the points x in the disjoint intervals 
(ui, v\), . . . , (uff, uat). An expression for the the reduced density matrix pa is obtained 
from (7) by sewing together only those points x which are not in A. This has the effect 
of leaving open cuts, one for each interval (v,j, Vj), along the line r — as in the right 
panel of figure 1. 

We may then compute Tr p n A , for any positive integer n, by making n copies of the 
above, labelled by an integer j with 1 < j < n, and sewing them together cyclically 
along the cuts so that <pj (x, r = f3~) = <fij+i(x, r — + ) and § n (x,r = (3~) = 4>\(x,t = 
+ ) for all x G A. This defines an n-sheeted structure depicted for n — 3 and in the 
case when A is a single interval in Fig. 2. The partition function on this surface will 
be denoted by Z n (A) and so 

^Pa = ^- (8) 

When the right hand side of the above equation has a unique analytic continuation to 
Ren > 1, its first derivative at n = 1 gives the required entropy 

S , = -h4^. (9) 
n^i dn Z n W 
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Figure 2. A representation of the Riemann surface 7^3,1. Reprinted with 
permission from [12]. 

So far, everything has been for a discrete space domain. We now discuss the 
continuum limit, in which a — ► keeping all other lengths fixed. The points x then 
assume real values, and the path integral is over fields <p{x 1 r) on an n-sheeted Riemann 
surface, with branch points at Uj and Vj. In this limit, Se is supposed to go over into 
the euclidean action for a quantum field theory. We indicate these n-sheeted surfaces 
with TZ n ,N and they are fully defined by the 2N branch points Uj and Vj. Whenever 
the value of n and N is not important, we will simply indicate the surface with 1Z. 

In the following, we will restrict our attention to the case when the quantum 
field theory is Lorentz invariant, since the full power of relativistic field theory can 
then be brought to bear. The behaviour of partition functions in this limit has been 
well studied. In two dimensions, the logarithm of a general partition function Z in a 
domain with total area A and with boundaries of total length C behaves as 

log Z = hAa- 2 + f 2 Ca- 1 + ... (10) 

where f\ and fa are the non-universal bulk and boundary free energies. Note, however, 
that these leading terms cancel in the ratio of partition functions in (8). 

In a conformal field theory, as was argued by Cardy and Peschel [14], there are also 
universal terms proportional to log a. These arise from points of non-zero curvature 
of the manifold and its boundary. In our case, these are conical singularities at the 
branch points. In fact, it is precisely these logarithmic terms which give rise to the 
non-trivial dependence of the final result for the entropy on the short-distance cut- 
off a. For the moment let us simply remark that, in order to achieve a finite limit 
as a — + 0, the right hand side of (8) should be multiplied by some renormalization 
constant Z(A,n). 

2-4- From replicated world-sheet to replicated target-space: Twist fields 

In the simplest instances it is possible to directly calculate the partition function on a 
n-sheeted Riemann surface, but in most of the cases this is very difficult. However, the 
surface we are dealing with has curvature zero everywhere except at a finite number 
of points (i.e. the boundaries between A and B Uj,Vj above). Since the lagrangian 
density does not depend explicitly on the Riemann surface 1Z as a consequence of 
its locality, it is expected that the partition function can be expressed as an object 
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calculated from a model on the complex plane C, where the structure of the Riemann 
surface is implemented through appropriate boundary conditions around the points 
with non-zero curvature. Consider for instance the simple Riemann surface lZ n _i 
needed for the calculation of the entanglement entropy of a single interval 
made of n sheets sequentially joined to each other on the segment x £ [tti, Vi], t = 0. 
We expect that the associated partition function in a theory defined on the complex 
plane z = x + it can be written in terms of certain "fields" at z = v\ and z = u\. 
The partition function (here £[tp](x, r) is the local lagrangian density) 



(11) 



(12) 



dxdr C[ip] (x, t) , 

m 

essentially defines these fields, i.e. it gives their correlation functions, up to a 
normalization independent of their positions. However in the model on the complex 
plane, this definition makes them non-local (see for a complete discussion [12]). 
Locality is at the basis of most of the results in field theory, so it is important to 
recover it. 

The solution to the problem consists in moving the complicated topology of the 
world-sheet 1Z (i.e. the space where the coordinates x, r live) to the target space (i.e. 
the space where the fields live). Let us consider a model formed by n independent 
copies of the original model. (Note that n is the number of Riemann sheets necessary 
to describe the Riemann surface by coordinates on the plane.) The partition function 
(11) can be re- written as the path integral on the complex plane 

Z n = / [dipi ■ ■ ■ dip„] exp - / dxdr (C[(fi](x,T) + . . . + C[ip n ](x,T)) 
Jc ultV1 L Jc 

where with \ r we indicated the restricted path integral with conditions 

tpi(x,0 + ) = Lp i+1 (x,0~) , x e i = l,...,n (13) 

where we identify n + i = i. The lagrangian density of the multi-copy model is 

£ (n) [<Pl,...,<P n ](x,T)=£[<p 1 ](x,T) + ... + £[(fi n } (x, t) 

so that the energy density is the sum of the energy densities of the n individual copies. 
Hence the expression (12) does indeed define local fields at (ui,0) and («i,0) in the 
multi-copy model [12]. 

The local fields defined in (12) are examples of "twist fields". Twist fields 
exist in a quantum field theory whenever there is a global internal symmetry a (a 
symmetry that acts the same way everywhere in space, and that does not change 
the positions of fields): J dxdr C[<jip](x,T) = J dxdr £[(p] (x, r). In the model with 
lagrangian there is a symmetry under exchange of the copies. The twist fields 

defined by (12), which have been called branch-point twist fields [12], are twist fields 
associated to the two opposite cyclic permutation symmetries i ^ i + 1 and j + l^i 
(i = 1, . . . , n, n + 1 = 1). We can denote them simply by T n and T n , respectively 

T n = T a , a : in i + 1 mod n , (14) 

7~ n = T a -i , cr _1 : i+lHimodn. (15) 

Notice that T n can be identified with 71„ (and in fact they were called $„ and $_„ 
in Ref. [3]). 

For the n-sheeted Riemann surface along the set A made of N disjoint intervals 
[iij,Vj] we then have 

Zn n , N cx (T n («i,0)t n («i,0)"-7;(« JV ,0)t n («j V ,0)) jC („) >c . (16) 
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Figure 3. Uniformizing transformation for Tt n i. ui — > £ = (ui — u)/(w — v) maps 
the branch points to (0, oo). This is uniformized by the mapping f — » 2 = £ 



This can be seen by observing that for x G consecutive copies are connected 

through t — due to the presence of T n (vj , 0), whereas for x in B, copies are connected 
to themselves through r = because the conditions arising from the definition of 
T n (uj,0) and 7~ n (vj,0) cancel each other. More generally, the identification holds for 
correlation functions in the model C on TZ n> \ 

lr ., i , ■ \ . (T n {u 1 ,Q)f n {v u Q)O l {x 1 T)---) c[nKC 

{0{x, r; sheet t)-. •) £ , 7 e n , 1 = — — p- — (17) 

where Oi is the field in the model coming from the i th copy of C, and the ratio 
properly takes into account all the proportionality constants. The same expression 
with the products of more twist and anti-twist fields holds in the case of TZn.u- 
It is also useful to introduce the linear combination of the basic fields 

n 

£ fc = £ e 2 "**^ , k = 0,1,..., n-1, (18) 

3 = 1 

that get multiplied by e 2lxlk l n on going around the twist operator, i.e. they diagonalizc 
the twist 

T n ^ k = e 2 " fc /'V fc , and f n <p k = e^/"^ . (19) 

Notice that when the basic field ipj are real (p* k — (p n -k- When the different values of 
k decouple, the total partition function is a product of the partition functions for each 
k. Thus also the twist fields can be written as products of fields acting only on (p k 

n — 1 n—1 

%=Y[ Tn,k , f n = J] Tn,k , (20) 

fc=0 fe=0 

with T fc! „<^ fe / = (p k iik^k' and Tk, n <Pk = e 2lx%k l n (p k . Thus 

n-1 

,n{vi,0) . . ■)£(»), C • ( 21 ) 

This way of proceeding is useful for free theories as in Refs. [15, 16, 17, 18, 19], when 
the various fc-modes decouple leading to Eq. (21). 



3. Entanglement entropy in conformal field theory: a single interval 

Following Ref. [3], we first consider the case N = 1 and no boundaries, that is the case 
considered by Holzhey et al. [1] of a single interval [u,v] of length I = \u — v\ in an 
infinitely long ID quantum system, at zero temperature. The complex coordinate is 
w = x + it and w = x — it. The conformal mapping w — > £ = (w — u)/(w — v) maps 
the branch points to (0,oo). This is uniformized by the mapping £ — > z — C}l n = 
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({w-u)/(w-v)) 1/n . This maps the whole of the n-sheeted Riemann surface 7£„ ; i to 
the z-plane C, see Fig. 3 for an illustration of this. Now consider the holomorphic 
component of the stress tensor T(w). This is related to the transformed stress tensor 
T(z) by [20] 

TM=(±L) 2 nz) + ±{zM, (22) 

where {z, w} — (z'"z ! — |z" 2 )/z' 2 is the Schwarzian derivative. In particular, taking 
the expectation value of this, and using (T(z))c — by translational and rotational 
invariance, we find 



<^-&».»}- *;;V ■)>-.)- ■ ,23) 

From Eq. (17), this is equal to 



(T n (u,0)f n (v,0)T J (w)) cinhc 
(T„(w,0)T„(w,0)) £ („) iC ' 

for all j. We can then obtain the correlation function involving the stress-energy tensor 
of by multiplying by n: 

(T n (u,0)f n (v,0)T^(w)) cln ^ c _ c (n 2 - 1) (u - v) 2 



(T n (u,0)T n (v,0)) cMtC 24n (w - u) 2 {w - v) 2 ' 

The comparison with the conformal Ward identity [20] 
(T n (u,0)f n (v,0)T^(w)) cMtC = 

1 d hr 1 d h-f \ 

+ 7 -^2+ a~ + , \ 2 ) (T n (u,0)T n {v,0)} c( n )tC , (24) 



w — udu (w — u) 2 w — v dv (w — v)' 

allows us to identify the scaling dimension of the primary fields T n and T n (they have 
the same scaling dimension d n = d n ) using (T n (u, 0)T n (v, 0)) £ (n) c = \u— v\~ 2drl | 

**~±(n-l). (25) 

To our knowledge this scaling dimension was first derived by Knizhnik [21] in a different 
context. 

The above equation determines all the properties under conformal transforma- 
tions, and we therefore conclude that the renormalized Z n (A)/Z n cx Trp^ behaves 
(apart from a possible overall constant) under scale and conformal transformations 
identically to the two-point function of a primary operator with dimension d n . In 
particular, this means that 

Trp" A = Cn ^—) . (26) 

The power of a (corresponding to the renormalization constant Z) has been inserted 
to make the final result dimensionless, as it should be. The constants c n cannot be 

$ We use here d n instead of A„ in [3] to avoid confusion between these scaling dimensions, in fact 
they are not the same and are related by d n = 2nA„. 
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determined with this method. However C\ must be unity. The analytic continuation 
is straightforward leading to the Renyi and von Neumann entropies 



where we defined the non universal constant 

, _ logC„ 



(28) 



n l-n 

Notice that c[ is minus the derivative of c n at n = 1. Despite their non- universal 
nature, the constants c' n are known exactly for few integrable models [22, 23, 24, 25, 12]. 

Under certain conditions, the entanglement entropy can also be expressed in terms 
of averages over ensembles of random matrices [26] providing a new connection between 
the universality class of the conformal field theory and random matrix ensembles. 

Interestingly, Eq. (26) describes the asymptotic behaviour for large enough £ for 
any n > (and not only for n > 1) at least for the simplest solvable models, where it 
can be explicitly checked. This is of relevance for the convergence of of the algorithms 
based on matrix product states [27, 28]. 

3.1. Generalizations: Finite temperature or finite size 

The fact that Trp^| transforms under a general conformal transformation as a two- 
point function of primary operators T means that it is simple to compute in other 
geometries, obtained by a conformal mapping z — ► w = w(z), using the formula 

(T n {z 1 ,z 1 )f n (z 2 ,z 2 )) = \w\z 1 )w\z 2 )\ dn (T n {w 1 ,w 1 )f n {w2,w 2 )}. (29) 

Particularly relevant is the transformation w — > z = {[3 /2ir) log w that maps each 
sheet in the w-plane into an infinitely long cylinder of circumference f3. The sheets 
are now sewn together along a branch cut joining the images of the points u and v. 
By arranging this to lie parallel to the axis of the cylinder, we get an expression for 
Tr p\ in a thermal mixed state at finite temperature (i~ x . After simple algebra, this 
leads to the result for the entropy [3, 29] 



Sa = \ log ( — sinh ^- ) + c\ = < 
3 \ira p 



c i 

- log - + c[ t<.0, 
6 a 

:£ + c[ 



(30) 



3/3 

This simple formula interpolates between the zero-temperature result for I <C fi and 
an extensive form in opposite limit I ^> (3. In this limit its density agrees with that 
of the Gibbs entropy of an isolated system of length £, as obtained from the standard 
CFT expression [30, 31] (3F <~ —(irc/6)(£/0) for its free energy. As expected, in the 
high-temperature limit, the von Neumann entropy reduces to a pure thermal form and 
the entanglement vanishes. 

Alternatively, we may orient the branch cut perpendicular to the axis of the 
cylinder, which, with the replacement (3 — > L, corresponds to the entropy of a 
subsystem of length I in a finite ID system of length L, with periodic boundary 
conditions, in its ground state. This gives 

fL. ^\- C( "~ 1/rl)/6 c (L . n£\ , 

Trp^ = c„ — sin— , S A = - log — sin — + c[ . (31) 

\ 7ra L I 3 \ ira L J 
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Figure 4. Exact fi- 
nite size scaling of the 
entanglement entropy 
S„(L) (here n = I) 
for the XXZ model 
at A = 1/2 against 
the CFT prediction 
1/3 log sin(-7rn / L) + 
(full line). c' x = 
0.7305 has been fixed 
[33]. Reprinted with 
permission from [32] . 



Sa is symmetric under t — > L — £ and it is maximal when £ = L/2. This relation 
is fundamental in the analysis of numerical data that are mainly done in finite size. 
It provides an unambiguous way to determine the central charge even from relatively 
small system sizes. In Fig. 4 we report the exact calculation of Sa for the XXZ chain at 
A = 1/2 from Ref. [32], showing that already small values of £ (< 6) it gives the correct 
asymptotic scaling. We stress that the most powerful aspect in the determination of 
the central charge via the entanglement entropy is that it does not involve the a priori 
knowledge of the speed of the sound, unlike other measures based on the gap or free 
energy scaling. For this reason, this method has been widely used in recent years. 



3.1.1. Finite temperature and finite size. When a finite system is also at finite 
temperature, we need to consider periodic boundary conditions both on space and 
imaginary time axes. This results in the calculation of a two-point function of twist 
operators on a torus. As it is well known [34], in this case correlations depend not 
only on the scaling dimensions, but on the full operator content of the theory and the 
calculations must be done for any universality class. Furthermore it is not possible to 
use a uniformizing transformation from the plane, because of the non-trivial topology 
of the torus. 

To our knowledge, this calculation has been performed only for the massless Dirac 
fcrmion [35]. In this case, it is convenient to use the representation of the twist fields 
in the diagonal basis (see Eq. (21)) to obtain 

n— 1 

Trp n A = ]J(T ntk (z,z)f ntk {0,0)) , (32) 

k=0 

and these two-point functions are known from bosonization (setting z = t and L = 1 
as scale) 



(T n>k {t)f n , k {0)) 



2irri(if3) :: 



2fc 2 /» 



9„(k£/n\i/3)\ 2 



(33) 



0i(£\i/3) 

where 6*„ and r\ represent standard elliptic functions, and v = 2, 3, 4 is the sector of the 
fermion. This gives a simple and compact answer for any integer n, but the analytic 
continuation of the second part is complicated because fc enters in the argument 
of the 9 functions. For this reason, it is possible to give exact expressions for the 
entanglement entropy only in the high- and low-temperature expansions (that however 
give convergent expressions). As an example we report here the high-temperature 
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expansion in the NS sector [y = 3) from Ref. [35] 



Sa = g log 



3 . n£ 
— sinn — 
na p 



1 ~ (l_ e 2^ e -2^ )(1 _ e -2^ e 



-2tt 



+ 5l> 



7 



3^ (l-e- 2 ^)2 



oo , y ZL« C oth ( 2# - 1 

;=i sinh (jrjj 

This formula gives one example of the crossover form Eq. (31) for » 1 to Eq. (30) 
for <C 1. More details about the derivation and the results for other sectors can be 
found in Ref. [35]. We mention that these issues have been investigated numerically 
for different spin-chains in Ref. [36]. 

3.2. Systems with boundaries 

In numerical simulations with DMRG and in real experimental life, physical systems do 
not obey periodic boundary conditions, but rather have some open boundaries. While 
in the study of correlation functions of local operators the effect of the boundaries can 
be reduced by performing measures far from them, the intrinsically global nature of 
the block entanglement makes it more sensible to the boundary conditions. This is 
not a negative feature and can be effectively described by boundary CFT [37, 38] . 

Let us start by considering a ID system is a semi-infinite line, say [0,oo), and 
the subsystem A is the finite interval [0,£). The n-sheeted Riemann surface TZ n ,i then 
consists of n copies of the half-plane x > 0, sewn together along < x < £, t = 0. 
Once again, we work initially at zero temperature. It is convenient to use the complex 
variable w = r+ix. The uniformizing transformation is now z = {{w— i£)/ (w+ii)*) 1 ^™ , 
which maps the whole Riemann surface to the unit disc \z\ < 1. In this geometry, 
(T(z)) — by rotational invariance, so that, using (22), we find 

'^.■4"-»"> (,-i)> ti ,f < 35) 

Note that in the half-plane, T and T are related by analytic continuation: T(w) = 
[T(w)]* [37]. Eq. (35) has the same form as (T(w)T n (i£)), which follows from the 
Ward identities of boundary CFT [37], with the normalization (T n (i£)) = (2£)~ dn . 
The analysis then proceeds in analogy with the previous case leading to 

2£\ ( c / 12 )("- 1 /«) c 2£ 



Trp n A = c n - => S A = -log-+c' 1 . (36) 

\ a J 6 a 

The constants c„ are in principle different from c„ in the periodic case. The coefficient 
in front of the logarithm is one-half of the one with periodic boundary conditions. This 
can be interpreted as the analogue of the area law in ID. In fact, while with periodic 
conditions there are two boundary-points between A and B, in the present case there 
is only one. 

Once again, this result can be conformally transformed into a number of other 
cases. At finite temperature we find 



c 



2£ 



( 8 2ir£\ I ^ log — + c\ £</3 , 

— £ + £[ £->d, 



6 \7ra (3 
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In the limit f >/3we find the same extensive entropy as before. This allows to identify 
[3, 40, 43] 

5i-y=log<7, (38) 

where log g is the boundary entropy, first discussed by Affleck and Ludwig [39]. g 
depends only on the boundary CFT and its value is known in the simplest cases. 
Numerical simulations confirm with high precision this relation [40]. It is worth 
mentioning that the change in the entanglement entropy of topological quantum 
Hall fluids (see the Fradkin review in this special issue [41] for details) at a point 
of constriction is related to the change of Affleck and Ludwig entropy of the coupled 
edge states of the fluid at the point contact [42] . 

For a finite ID system, of length L, at zero temperature, divided into two pieces 
of lengths t and L — t, we similarly find 

S A = - log — sin — + c[ . 39 
6 \7ra L J 

This last equation is the most appropriate for numerical simulations that are usually 
performed in finite systems with some boundary conditions at both ends. 



3.2.1. Interfaces. We have seen that when a system is translationally invariant the 
entanglement entropy scales like Sa = c/31og^, while in the presence of a boundary, 
that can be a disconnected chain, it scales like Sa = c/61og^. In a condensed matter 
system like a spin chain, we can modulate a single bond (let say at x = 0) from zero 
to the value in the rest of the chain, going from one extreme to the other. In the 
presence of such a defect, there are mainly three possibilities under renormalization 
group: 

• The defect is irrelevant: the system flows to the translational invariant 
Hamiltonian and Sa = c/31og£. 

• The defect is relevant: the RG flow leads the system to a different fixed point. In 
particular when non-trivial ones are excluded, it flows to a disconnected system 
with S A = c/61og£. 

• The defect is marginal: the critical properties, and in particular the entanglement 
entropy, are continuous function of the defect strength. 

It has been shown numerically [44] and analytically [45] that in the gapless phase 
of the XXZ chain with A ^ 0, the defect is either relevant or irrelevant, leading always 
to the well-known behaviors of Sa- More interesting is the case of the XX model [46], 
when the defect is marginal and for the entanglement entropy one gets 

S A = ^log£, (40) 

where t is the strength of the defect (t = for disconnected chains and t = 1 
for translational invariant ones), and o~(t) is a monotonous increasing function of t 
with ct(0) = 1/2 and cr(l) = 1. A similar behaviour has been found also for more 
complicated defects in Ref. [47]. 

This phenomenon can be described as an interface between two different CFT's 
with c = 1. In Ref. [48] the entanglement entropy of two systems of length L separated 
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by an interface with scattering amplitude s (the analogous of t above) has been 
calculated 



(the integral can also be written in terms of polylog functions). When there is no 
interface, i.e. for s — 1, cr(l) = 1/3, as expected. Instead cr(0) = because the 
two CFT's decouple. Unfortunately no result for a subsystem of length I < L in the 
presence of the interface is still available to be compared with the results in the XX 
chain [46]. 

Other results for more general defects are known [49], but we remand to the 
review by Laflorencie et al. in this issue [50] for an extensive discussion. 

3. 3. General appearance of logarithmic behaviour 

In arbitrary dimension, the scale invariance (i.e. criticality of the statistical model) 
together with translational and rotational (i.e. Lorentz in real time) invariance and 
locality automatically leads to conformal invariance [34], explaining the very large 
interest in these theories. However, nature is not always so kind and there are physical 
systems that are critical, but that because of the explicit breaking of translational 
and/or rotational invariance are not conformal. It is then natural to ask what is the 
behaviour of the entanglement entropy in these systems. Srednicki [9] argued that the 
area law in higher dimensional systems for a gapless model should generally collapse 
to a log^ behaviour in one dimension, and so one would expect the appearance of 
logarithms even in non-conformal invariant systems. Unfortunately, nowadays there 
are several examples of the breaking of the area scaling in critical systems (see e.g. 
Ref. [6]), leaving doubts on the earlier argument. 

Critical systems showing the breaking of translational or rotational invariance 
have been largely studied. We give here few important examples. When translational 
invariance is broken by quenched disorder, the resulting statistical model can be 
studied by strong-disorder RG methods and numerically. In all the studied models, it 
has been shown unambiguously that the entanglement entropy always shows a log^ 
behaviour [51, 52, 53, 54, 55, 56, 57, 58, 59] (see the review by Moore and Refael [60] in 
this issue). Translational invariance can also be broken by taking aperiodic couplings: 
even in this case a logt behaviour has been found [61, 62]. Non-relativistic dispersion 
relations like E = k 2 also naturally breaks conformal invariance, by breaking Lorentz. 
A well-known and physical important example is the ferromagnetic Heisenberg chain, 
for which the entanglement entropy scales like Sa — 1/2 log £ [63, 64]. Another 
interesting example of this kind can be found in [65]. 

Often it has been proposed that the scaling of the entanglement entropy as log£ 
can be used to define an effective central charge for non conformally invariant systems 
that can share some of the properties of the real central charge (as for example the 
monotonicity along renormalization group flow [66]). However, conformal invariance 
is so powerful that fixes the scaling of the entanglement entropy, but, as we have 
seen, it also gives precise predictions for finite size scaling in Eq. (31) and for the 
scaling of the Rcnyi entropies (26). Before arguing about asymptotic restoration of 
conformal invariance, all these relations should be carefully checked. For example the 
finite size scaling found in Ref. [64, 65] is different from Eq. (31). Also the scaling 
of the entanglement entropy in the zero temperature mixed state of the XXZ chain 
at A = 1/2 scales logarithmically, but has a finite size form different from Eq. (31) 




(41) 
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[32] and cannot be described by CFT. Some excited states in spin-chains also displays 
logarithmic behavior because they can be interpreted as ground-states of properly 
defined conformal Hamiltonian [67]. For random systems with quenched disorder, 
the finite size scaling form seems to be conformal from numerical simulations, but 
the Renyi entropies have different scaling in the random singlet phase [60]. We also 
mention that in several collective models a similar logarithmic behavior has been found 
[68] , but its origin is different from the one discussed here because of the absence of a 
spatial structure in these mean-field like models. Their properties are related to those 
of the particle partitioning reviewed by Haque et al. in this volume [69] . 



4. Entanglement of disjoint intervals 

When the subsystem A consists of several disjoint intervals, the analysis becomes 
more complicated. In Rcf. [3] we provided a result that in general is incorrect. This 
was based on a uniformizing transformation mapping !Z n .N into the complex plane. 
However, the surface 1Z n ^N has genus (n — l)(N — 1) and so for N ^ 1 (that is the 
case we already discussed) cannot be uniformized to the complex plane (at the level of 
the transformation itself, this has been discussed in details [70] ) . The case n = N = 2 
has the topology of a torus, whose partition function depends on the whole operator 
content of the theory and not only on the central charge. Consequently the simple 
formulas of Ref. [3] cannot be generally correct. The partition functions on Ricmann 
surfaces with higher genus are even more complicated. 

Let us then start our analysis from the case of two intervals given by the surface 
By global conformal invariance the partition function (that is the four-point 
correlation of twist fields) can be written in the form (in this section we adsorb the 
normalization Z n into Zn n N ) 

Tr P \ = Z Un , 2 = d (-. \^-^-V2\ \ M " 

- vi\\u 2 - v 2 \\u 1 - v 2 \\u 2 - vi\J 

where x is the real four-point ratio 

= Z12 zu _ (ui - vi){u 2 - v 2 ) 
_ 213^24 («1 - U 2 ){V 1 - V 2 ) ' 

and d n is given by Eq. (25). It can also be written as 

v 2 =^/-w> ( 44 ) 

where Z^ is the incorrect result in [3]. We normalized such that T n (Q) = 1 (for 
x — > 0, Zfi n 2 is the product of the two two-point functions previously calculated 
and normalized with c n ). The function F n (x) depends explicitly on the full operator 
content of the theory and must be calculated case by case. 

In Ref. [71], using old results of CFT on orbifolded space [72, 73], T 2 {x) has been 
calculated for the Luttinger liquid CFT, that is a free boson compactified on a circle 
of radius R 

- Mt)]2 > ( 45 ) 

where r is pure-imaginary, and is related to x via x = [6*2(t)/^3(t)] 4 - v are Jacobi 
theta functions, rj is proportional to the square of the compactification radius R (while 
the definition of R usually depends on the normalization of the action, 77 is the same in 
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Figure 5. The mutual information for fixed four-point ratio x as function of r\ 
in the gapless phase of the XXZ model. The horizontal lines stand for ^ of 
Ref. [3]. Left: mutual information of the von Neumann entropy. Right: mutual 
information of the Renyi entropy for n = 2, compared with the compactified boson 
prediction Eq. (45). Reprinted with permission from [71] 



all literature, that is why we prefer to write everything only in terms of r\ that usually 

is written as r\ — l/(2if) in Luttinger liquid notation). 

In order to check this prediction and the failure of , in Ref. [71] the 

entanglement of the XXZ chain for generic values of the anisotropy A and magnetic 

field always in the gapless phase has been calculated by direct diagonalization for 

systems up to 30 spins. In the absence of the magnetic field, ij is related to the 

anisotropy by r\ = 1 — (arccos A)/7r, while for non-zero h z a closed formula for r\ does 

not exist and must be calculated numerically as explained in [71]. The main results 

coming from the exact diagonalization are reported in Fig. 5 in terms of the Renyi 

mutual information 

r(") _ c(") I c(") _ o(«) (AR\ 
A\.Ai Ai + D A 2 D AiUA 2 ' 

where A± and Ai are the two intervals composing A = A\ U Ai. In Fig. 5 the mutual 
information for n = 1,2 is reported. The r\ dependence in both the cases is evident, 
and the good collapse of the data shows the correctness of the scaling in Eq. (42). 

(2) 

(In I AvA ^ the collapse is worst due to the oscillating corrections to the scaling of 
the Renyi entropies, known already for the single interval [33]). In the right panel of 
Fig. 5, the comparison of the numerical results with the prediction (45) is reported. 
The agreement is satisfactory, considering the small subsystem sizes and the strong 
oscillations. The results for i^ .j^ in the left panel of Fig. 5 are more stable because 
of the absence of oscillations. These are then more suitable for the comparison with 
an analytic calculation. 

In Ref. [74] (in collaboration with E. Tonni), we managed to calculate T n {x) for 
generic integral n > 1. The result reads 

^W ^' ""'' 9 ^ . (47) 

e(o|r) 

where T is an (n — 1) x [n — 1) matrix with elements 

(48) 



2i v— \ / k , 
T rs = — sin [ 7r- j fj k/n cos 



k = l 



n 



2n—{r — s) 



n 
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and 

& = Fv vm X) ' ^ 2^1(2/, i-y;i;^) • (49) 

rj is exactly the same as above, while 9 is the Riemann-Siegel theta function 

0(z|r) = ^ exp [iirm ■ T ■ m + 2mm ■ z] , (50) 

that for n — 1 = 1 reduces to the Jacobi 6* 3 (r = i/?i/ 2 ), reproducing Eq. (45). 

Unfortunately we have been not yet able to continue analytically this result to 
real n for generic 77 and so obtain the entanglement entropy to compare with the left 
panel of Fig. 5 from Ref. [71]. 

Some interesting properties can be readily deduced from Eq. (47): 

• For any n, it is symmetric under the exchange 77 «-> I/77, generalizing the result 
for n = 2 in Eq. (45). 

• For 77 = 1 the result Z^ n 2 is correct, i.e. T n {x) = 1. 

• It is symmetric under the exchange x <-> 1 — x (crossing symmetry). This the 
scaling limit of Sa = Sb for finite systems [71]. 

These three important properties carry over to the analytic continuation and so must 
be true for the Von-Neumann entropy at n = 1. These findings then explain the 
numerical results of Ref. [71] (in the left panel of Fig. 5 the symmetry r\ <-> I/77 and 
T' x {x) = 1 are evident). 

More details and other interesting properties can be found in Ref. [74] . We discuss 
here the so-called uncompactified limit for rj » 1 or by symmetry tj<l. In this case 
we have [74] for 77 <C 1 

v -(n-l)/2 v -(n-l)/2 

Tn[X) ^ [Y\UF k/n {x)F k/n{ l-x)f^ [D n (x)D n (l-x W /z> (51) 

where the function D n (x) — Hfe=i Fk/n( x ) nas been analytically continued in Ref. [74]. 
Then the prediction 

C. A2 (v « i) - i%;Z - 4 ^77 + EM±Wz*l , (52) 

perfectly agrees with the numerical results in Fig. 5 (from Ref. [71]) for 77 < 0.4 (again 
in the previous equation _Y is the result of Ref. [3] ) . Also the regime i<1 can 
be studied analytically [74]. 

A few comments are now in order. In Ref. [16, 75, 76], the entanglement 
entropy for two disjoint intervals has been calculated for a free the Dirac fermion, 
that corresponds to a compactified boson with 77 = 1/2 [34]. However, it has been 
found that the entanglement entropy is given by , in contrast with the numerical 
calculations in Ref. [71] and the analytic one in [74]. The details of this apparent 
disagreement are still not completely understood, but they should be traced back to 
the different boundary conditions that result from constructing the reduced density 
matrix for spin or fermion variables. Another calculation in agreement with 2 can 
be found in [77]. For the Ising model numerical computations [78] also show a good 
agreement with . Also in this case, it is likely that the deviations from Z^ 

should be attributed to the choice of the variables used in constructing the reduced 
density matrix. (In fact, calculations in the spin variables [79] show numerically and 
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analytically that Z^ is not correct.) Finally, holographic calculations in AdS/CFT 
correspondence [80, 81], considering the classical limit in the gravity sector, also found 

2 . It would be interesting to understand how the correct result might arise from 
taking into account the quantum effects on the gravity side (for details see the review 
of Nishioka, Ryu, and Takayanagi in this volume [82]). 

For general N > 2, there are still no firm results in the literature. By global 
conformal invariance one can deduce 

/rr ( \( \ \ ( c / 6 )(™- 1 /") 

**-*( 't,^ ) <53) 

For .Fn,jv({a;}) = 1 this is the incorrect result of Ref. [3] (note a typo in the 
denominator), {x} stands for the collection of 2N — 3 independent ratios that can 
be built with 2N points. Some old results from CFT on orbifold in Refs. [73, 83] 
could be useful to calculate !Fn,N{{x}) for a compactified boson. We also mention 
that in 2D systems with conformal invariant wave-function (reviewed in [41] in this 
volume), the entanglement entropy of a single region displays an addictive universal 
term depending on 77 [84]. 

Finally, it is worth to recall that in the case of more intervals, the entanglement 
entropy measures only the entanglement of the intervals with the rest of the system. 
It is not a measure of the entanglement of one interval with respect to the others, that 
instead requires the introduction of more complicated quantities because A\ U A2 is 
in a mixed state (see e.g. Refs. [85] for a discussion of this and examples). 



5. Entanglement entropy in non-critical 1+1-dimensional models 

When a one-dimensional statistical model has a gap (i.e. the underlying quantum field 
theory is massive) the entanglement entropy saturates to a finite value [2]. This is an 
analogue of the area law in one dimension, because the area is only a number that 
does not increase with subsystem size, in contrast to higher dimensions. Generally 
this value is very complicated to be calculated and it is known only in very simple 
cases. However, when a system is close to a conformal quantum critical point, that is 
when the gap A is small (or the correlation length ^ a A -1 is large) it is possible to 
derive a very general scaling form [3] , that can be used also as an operative definition 
of the correlation length. Hastings [86] (see also [87]) provided a rigorous proof of the 
area-law for one-dimensional systems with a generic local Hamiltonian, not necessarily 
close to a conformal critical point. 

We consider an infinite non-critical model in 1+1-dimcnsions, in the scaling limit 
where the lattice spacing a — » with the correlation length (inverse mass) fixed. This 
corresponds to a massive rclativistic QFT. We first consider the case when the subset 
A is the negative real axis, so that the appropriate Riemann surface has branch points 
of order n at and infinity. However, for the non-critical case, the branch point at 
infinity is unimportant: we should arrive at the same expression by considering a finite 
system whose length L is much greater than £. 

Our argument parallels that of Zamolodchikov [66] for the proof of his famous 
c-theorem. Let us consider the expectation value of the stress tensor of a massive 
euclidean QFT on such a Riemann surface. In complex coordinates, there are three 
non-zero components: T = T zz , T = T zz , and the trace 6 = 4T ZZ — 4T ZZ . These are 
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related by the conservation equations 

d g T+^d z 6 = and d z T+^d E Q = 0. (54) 

Consider the expectation values of these components. In the single-sheeted geometry, 
(T) and (T) both vanish, but (6) is constant and non- vanishing: it measures the 
explicit breaking of scale invariance in the non-critical system. In the n-sheetcd 
geometry, however, they all acquire a non-trivial spatial dependence. By rotational 
invariance about the origin, they have the form 

(T(z,z)) = F n (zz)/z 2 , (55) 

{Q(z,z))-(Q) 1 = G n {zz)/(zz), (56) 

(T{z,z)) = F n (zz)/z 2 . (57) 

From the conservation conditions (54) we have 



(zz) (f^ + \g>^ =\G n . 



(58) 



Now we expect that F n and G n both approach zero exponentially fast for \z\ >• £, 
while in the opposite limit, on distance scales <C £, they approach the CFT values 
F„^(c/24)(l-n- 2 ), G„^0. 

We define an effective C-function 

C n {R 2 )^{F n {R 2 ) + - A G n (R 2 )^ R 2 -JLj Cn (R 2 ) = ±G n (R 2 ). (59) 

whose integral, assuming that theory is trivial in the infrared (if the RG flow is towards 
a non-trivial theory, c should be replaced by cuv — cm), gives 

G n{R 2 ) M v2,__ C/ 1 



or cquivalcntly 



s: 
i 



r >2 d(R 2 ) = --[l-^) , (60) 



((e)n-(Q)i)d 2 R = -nn^(l-±^j , (61) 

where the integral is over the whole of the n-sheeted surface. Now this integral 
(multiplied by a factor 1/2-k corresponding to the conventional normalization of 
the stress tensor) measures the response of the free energy — log 2 to a scale 
transformation, i.e. to a change in the mass m, since this is the only dimensionful 
parameter of the renormalized theory. Thus the left hand side is equal to 

-(2^)mA[iogZ„-nlogZ] , (62) 

giving finally 

^ = cMW^M/n) ) (63) 

where c„ is a constant (with however c\ = 1), and we have inserted a power of a, 
corresponding to the renormalization constant Z discussed earlier, to make the result 
dimensionless. Differentiating at n = 1, we find 

c c £ 

Sa = - ~ log(ma) = - log - , (64) 
6 6 a 

where £ is the correlation length. We re-emphasize that this result was obtained only 
for the scaling limit £ 3> a. 



Entanglement entropy and CFT 



20 



So far we have considered the simplest geometry in the which set A and its 
complement B are semi-infinite intervals. The more general case, when A is a union 
of disjoint intervals, is more difficult in the massive case. However it is still true that 
the entropy can be expressed in terms of the derivative at n = 1 of correlators of twist 
operators T, T. The above calculation can be thought of in terms of the one-point 
function (T n ). In any quantum field theory a more general correlator (n»=i ®( w i)}: 
with <!> = T n or T n , should obey cluster decomposition: that is, for separations |iu» — Wj \ 
all 3> £, it should approach (T n } k . This suggests that, in this limit, the entropy should 
behave as 

S A = A%g^, (65) 
o a 

where A = k is the number of boundary points between A and its complement. This 
would be the ID version of the area law. When the interval lengths are of the order 
of £, we expect to see a complicated but universal scaling form for the cross-over. 

This scaling has been confirmed in several cases with A = 1 or 2 (see e.g. 
[3, 88, 23, 89, 90, 24, 91, 92, 93], but this list is far from being exhaustive). The 
corrections to this formula for I <C £ are also universal [12, 94, 95, 96] and are discussed 
in details in the review by Castro- Alvaredo and Doyon [97] in this issue. 



6. Entanglement spectrum 

The knowledge of the scaling form for Trp^ as in Eqs. (26,31,36,63) gives more 
information about the reduced density matrix than the entanglement entropy. We 
have seen that in many cases it scales like 

R n = Trp n A = c n e- b ^- 1 / n \ (66) 

with b > only depending on the main features of the set A, on the characteristic 
length of the system L e s (i.e. £, £, Lsmir£/L . . .) and on the central charge. This 
suggests that many properties of the reduced density matrix are very universal 
and do not depend on the details of the theory. For example, the scaling of the 
largest eigenvalues A max of pa is obtained by taking the limit for n — ► oo [98]: 
S A °°^ = — lnA max = Sa/2 defines the so-called "single-copy entanglement" [99] and 
gives another measure of the entanglement content of an extended system. This 
peculiar scaling led to the conclusion that "half the entanglement in critical systems 
is distillable from a single specimen" [100]. 

This result is however only the tip of a lot of information about the reduced 
density matrix obtainable from Eq. (66). This information is encoded in the full 
spectrum of the reduced density matrix, which has been shortly called "entanglement 
spectrum" [101] and has been derived in Ref. [102] for ID systems from Eq. (66). 

In order to characterize the entanglement spectrum, let us define the distribution 
of eigenvalues P(X) = J2 i S(X — Aj). If we ignore the n dependence of the coefficient c„ 
(that however is expected only to give corrections the leading behaviour) , it is easy to 
compute the entanglement spectrum by inverse Laplace transforming R n , obtaining 
[102] 

b0(X max — A) 



P(A) = S(X max - A) + - ^ / 7 1 (2 v /bln(A max /A)) , (67) 
A v /61n(A max /A) 

where Ik(x) stands for the modified Bessel function of the first kind. Amazingly, P(A) 
depends only on A max (we recall b = — lnA max ) and not on any other detail of the 
theory. For some integrable gapped systems P(A -C 1) was previously derived [103]. 




Figure 6. Sum of the first M eigenvalues of the XX model up to M = 100: 
1 - s(M) as function of M for t = 10, 100, 1000, 10000 (black dots). The red line 
is the conformal field theory prediction Eq. (70), in which A max has been fixed to 
the maximum eigenvalue obtained numerically. Numerical data from Ref. [102]. 



Let us discuss now the main properties of -P(A): 

The mean number of eigenvalues larger than a given A is 

n(A) = / ^ dXP(X) = 7 (2v / feln(A max /A)) . (68) 



Note that for A — » 0, n(A) diverges, as it should, because in the continuum the 
number of eigenvalues is infinite. In the lattice models, this can be regularized 
by the finite number of degrees of freedom. 

The normalization ^2 Ai = 1 corresponds to J \P(X)d\ = 1. 

The entanglement entropy is given by 

S = - / A In XP(X)d\ = -2 In A max , (69) 
Jo 

reproducing the result that the single copy entanglement equals one-half of the 
entanglement entropy. 

Majorization is a relation between two probability distributions A = {A.;} and 
/i = {/Xj} whose elements are ordered Ai > A2 • • • > Xn (and similarly for 
fi): it is said that A majorizes fi if J2i=i Ai > Si=i Mi f° r an Y M = 1, . . . , N 
and X^j=i A? = 2i=iA*j = ^ has been argued, observed numerically and 
in some instances proven analytically, that with increasing L e ff the resulting 
distribution of eigenvalues is majorized by the ones at smaller scaling lengths 
[104, 40]. Majorization follows easily from Eq. (67): 



M 

s{M) = ^2 Ai -> A max 

i=l 



L (M) 



1+/ dye" 2 ' /^(y) 



(70) 



at fixed M, is a monotonous function of A max (that is a monotonous function of 
L e ff)- This proves majorization. 

The ratio of the first two eigenvalues (see [102] for details) is 

■r—— = e k / b = e^^W^T ; (71) 
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where A: is a constant and in the second equality we used the result for periodic 
boundary conditions b — — In A max = (c/6) \n(£/a). This agrees with an old result 
for the scaling of eigenvalues of the corner transfer matrix [105]. 

Fig. 6 shows an explicit check of the distribution of eigenvalues for the XX model 
obtained in Ref. [102] by using the methods of Refs. [106, 107]. It is evident that 
when the subsystem size t is large enough, the numerical results perfectly agree with 
Eq. (70). 

The knowledge of the scaling form of the entanglement spectrum has been 
fundamental in understanding the convergence and the scaling of the algorithms 
based on matrix product states (MPS) [27], like DMRG. In these algorithms, the 
maximum amount of entanglement that can be effectively described is limited by 
the dimension \ of the matrix used to describe the state. The maximum possible 
entanglement of this state is 5 max = logx> when all the components have the same 
weight 1/x- But this maximum-entanglement state has nothing to do with the ground- 
state of the local Hamiltonian the algorithm is searching for, because it is described 
by Eq. (67). Numerical studies [108, 109], in fact, showed that the entanglement of 
the MPS approximating the critical ground-state scales like 

5 A = ylog X , (72) 

defining an effective length £ x ~ \ K [108]. k has been introduced as a new critical 
exponent of the MPS [108]. Using Eq. (68), it has been possible to calculate this 
exponent exactly [109], obtaining 

k = , 6 , (73) 

<yi2/ c +i - 

and to show that the corrections to Eq. (72) scale like l/log%. 

It is worth mentioning that a new numerical algorithm specifically based on 
the scaling properties of the entanglement in a conformal critical point has been 
recently proposed by Vidal [110]: the multi-scale entanglement renormalization ansatz 
(MERA). In a MERA, the ground-state of an extended quantum system is organized 
in layers corresponding to different length scales and, at a quantum critical point, each 
layer equally contributes to the entanglement of a block. This method then allowed 
to simulate systems of remarkably large sizes with a relatively little numerical effort. 
In particular, since the method is explicitly designed for scale invariant systems, some 
deep connections with CFT have been revealed [111]. 



7. Entanglement entropy after a quantum quench 

The experimental realization of cold atomic systems with a high degree of tunability of 
hamiltonian parameters, and the ability to evolve in time with negligible dissipation, is 
motivating an intensive study of extended quantum systems out of equilibrium. New 
numerical algorithms have been developed to describe the time-evolution of quantum 
systems effectively (among which adaptive time-dependent DMRG [112], known as 
tDMRG has been by now the most successful). As for the equilibrium counterpart, 
the amount of entanglement of the time-dependent state governs the effectiveness of 
the numerical methods based on MPS. It is then fundamental to have firm bases and 
expectations for the entanglement growth in non-equilibrium situations. In this case 
also conformal field theory has been a fundamental tool in understanding very general 
properties of the time evolution of the entanglement entropy. 
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Figure 7. Left: Space-imaginary time regions for the density matrix in (75). 
Right: The reduced density matrix pa is obtained by sewing together along r = 
only those parts of the x-axis corresponding to points in B (right part in this 
plot). 



The most studied situation (both theoretically and experimentally) is a so called 
quantum quench. In a quench, an extended quantum system is prepared at time t = 
in a pure state \ipo) which is the ground state of some hamiltonian H . For times 
t > the system is allowed to evolve unitarily according to the dynamics given by a 
different hamiltonian H , which may be related to H by varying a parameter such as 
an external field. 

Based mainly on the first results from conformal field theory [113, 114] it has been 
possible to understand that the entanglement entropy grows linearly with time for a 
so called global quench (i.e. when the initial state differs globally from the ground 
state and the excess of energy is extensive), while at most logarithmically for a local 
one (i.e. when the initial state has only a local difference with the ground state and 
so a small excess of energy) . As a consequence a local quench can easily be simulated 
by means of tDMRG, while a global one is harder and the numerics must be limited 
to relatively small systems sizes. New time-dependent algorithms based explicitly on 
the possibility of "storing" more entanglement [115, 116, 117] are being developed to 
have full access to these dynamics. 

7.1. Global quench 

Suppose we prepare the system in a state \ipo(x)) and unitarily evolve it with the 
hamiltonian H. The matrix elements of the density matrix at time t are 

(r{x")W)W{x')) = (r(xn\e- ltH \Mx))(Mx)\e +ltH W(x')) • (74) 
We modify this time dependent expectation value as 

{r{x")W)W{x')) = Z- 1 ^"(x")|e- ltff -^|^o(x))(^o(x)|e+^-^|^(x')) , (75) 

where we have included damping factors e~ eH in such a way as to make the path 
integral absolutely convergent. We shall see at the end of the calculation whether 
it is justified to remove them. The normalization factor Z = {ip (x)\e~ 2<iH \ipo{x)) 
guarantees that Tr p = 1 . 

Each of the factors may be represented by an analytically continued path integral 
in imaginary time: the first one over fields ip(x, t) which take the boundary values 
tf>o(x) on t = — 6 — it and ip"(x) on r = 0, and the second with i[)(z,t) taking the 
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values ip'{x) onr = and tpo( x ) on r = e — zt. This is illustrated in Fig. 7. Z is given 
by the euclidean path integral over imaginary time 2e, with initial and final conditions 
both equal to ipo(x). This is the same as sewing together the two edges in Fig. 7 along 
t = 0. As before, the reduced density matrix pA(t) is obtained by sewing together 
along t = only those parts of the x-axis corresponding to points in B, leaving open 
slits along A, and Tr p\ is given by sewing together n copies of this in a cyclic fashion. 
Thus 

Tr p n A = Z n {A)/Z n , (76) 

where Z n is the path integral on an n-sheeted surface, where the edges of each sheet 
correspond to imaginary times — T\ and t 2 , analytically continued to T\ = e + it and 
r 2 = e — it, and the branch points lie along r = at the boundaries points between 
A and B. Finally the entanglement entropy is given by the derivative of Trp^ with 
respect to n at n = 1. 

Eq. (75) has the form of the equilibrium expectation value in a 1 + 1-dimcnsional 
strip geometry with particular boundary conditions. We wish to study this in the 
limit when t and £ are much larger than the microscopic length and time scales, 
when RG theory can be applied. If H is at (or close to) a quantum critical point, 
the bulk properties of the critical theory are described by a bulk RG fixed point (or 
some relevant perturbation thereof). In that case, the boundary conditions flow to 
one of a number of possible boundary fixed points [118]. Thus, for the purpose of 
extracting the asymptotic behaviour, we may replace |^ ) by the appropriate RG- 
invariant boundary state \iPq) to which it flows. The difference may be taken into 
account, to leading order, by assuming that the RG-invariant boundary conditions 
arc not imposed at t = and r = 2e but at r = — t and r = 2e + t . In the 
language of boundary critical behaviour, To is called the extrapolation length [118]. 
It characterizes the RG distance of the actual boundary state from the RG-invariant 
one. It is always necessary because scale-invariant boundary states are not in fact 
normalizablc [38] . It is expected to be of the order of the typical time-scale of the 
dynamics near the ground state of Hq, that is the inverse gap m^ 1 . The effect of 
introducing To is simply to replace e by e + r . The limit e — > + can now safely be 
taken, so the width of the strip is then taken to be 2r . 



7.1.1. One interval in the infinite chain. Now we consider the case when H is critical 
and the field theory is a CFT. First let us consider the case when A is a slit of length 
£ and B is the rest of the real axis. For real r the strip geometry described above may 
be obtained from the upper half-plane by the conformal mapping w — (2r /n) log z, 
with the images of the slits lying along argz = 7tti/2to. The result for Z n /Z n in the 
upper half z-plane, with two branch points follows from Eq. (42) where half of the 
points are the images with respect to the real axis. To obtain the result in the strip 
geometry we transform this two-point correlation function according to Eq. (29). 

After some algebra and continuing to n = r +it (see [113, 119, 120] for detailed 
calculations), we find for t, £ » to 

LTPA(t) - Cn {2ro) \(e^o- e -^T o) 2 cosh 2 {nt/2To) ) -M*) • (77) 

T n (x) is the boundary analogue of T n (x) for the four-point function of twist fields in 
the plane in Eq. (42). After the conformal mapping and analytically continuing the 
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four-point ratio x becomes for £/t and t/r a large [119, 120] 

x jj^ -j— , (78) 

Thus for t < £/2 we have x ~ and in the opposite case t > 1/2 we have x <~ 1. 
Even if F n {x) is generally unknown, we only need its behaviour close to x <~ and 1, 
that are easily deduced from general scaling. Indeed when a; <~ 1 the two points are 
deep in the bulk, meaning !F n (l) = 1. Oppositely when the points are close to 

the boundary and again F n (0) = 1 (this because (T) ^ 0, see [120, 121]). Thus for 
the purpose of extracting the asymptotic behaviour, the function T n {x) is irrelevant, 
explaining why the results of Ref. [113], obtained within this assumption are correct. 
Putting everything together, in the case where £/tq and t/ro are large, the moments 
of the reduced density matrix simplifies to 

Trpm^c n (n/2r )^( ^±__ . (79) 



Differentiating wrt n to get the entropy, 

iret 



S A (t) ~ -^logr + 



6t 



t<£/2, 
t>£/2, 



(80) 



12r 

that is Sa (t) increases linearly until it saturates at t = IJ2. The sharp cusp in 
this asymptotic result is rounded over a region \t — £/2\ ^ t . As a difference with 
Ref. [113], following [122], we have added explicitly the subleading constant term log To 
confirming that To is connected to the inverse mass gap in the initial state. 

The result of the entanglement entropy for large time is the same of a mixed 
state at inverse large finite temperature (3 c s = 4to (see Eq. (30)). The physical 
interpretation of this important effect is that any finite subsystem A reaches a quasi- 
stationary thermal state, in which the infinite remaining part of the system B act as 
a thermal bath. It has been shown that, within CFT, this effective temperature is 
the same for any observable [120] and so can be properly defined. The possibility of 
defining a Gibbs-like asymptotic state for a general hamiltonian governing the time 
evolution (i.e. beyond the CFT case) is a subject of an intensive current activity that 
would require its own review and that will not be considered at all here. 

These results for translationally invariant states have been generalized to 
inhomogeneous quantum quenches with sharp [123] and smooth [122] initial states. 



7.1.2. Physical interpretation. The qualitative, and many of the quantitative, 
features of Sa{£) found above may be understood physically as follows [113]. The 
initial state |V>o) has a very high energy relative to the ground state of the hamiltonian 
H which governs the subsequent time evolution, and therefore acts as a source of 
quasiparticle excitations. Particles emitted from different points (further apart than 
the correlation length in the initial state oc To) are incoherent, but pairs of particles 
moving to the left or right from a given point are highly entangled. We suppose that 
the cross-section for producing such a pair of particles of momenta (p',p") is f(p',p"), 
and that, once they separate, they move classically. This will of course depends on H 
and the state I'i/'o), and in principle is calculable, but we made no strong assumptions 
on its form. If the quasiparticle dispersion relation is E = E pi the classical velocity is 
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2t>l 



Figure 8. Space- 
time picture illustrat- 
ing how the entangle- 
ment between an in- 
terval A and the rest 
of the system, due 
to oppositely mov- 
ing coherent quasipar- 
ticlcs, grows linearly 
and then saturates. 
The case where the 
particles move only 
along the light cones 
is shown here for clar- 
ity. Reprinted with 
permission from [113]. 



v p = dEp/dp. We assume that there is a maximum allowed speed which is taken to 
be 1, that is \v p \ < 1. A quasiparticle of momentum p produced at x is therefore at 
x + v p t at time t, ignoring scattering effects. 

Now consider these quasiparticles as they reach either A or B at time t. The 
field at some point x' e A will be entangled with that at a point x" € B if a pair 
of entangled particles emitted from a point x arrive simultaneously at x' and x" (see 
Fig. 8). The entanglement entropy between x' and x" is proportional to the length of 
the interval in x for which this can be satisfied. Thus the total entanglement entropy 
is 

S A {t)~ I dx' j dx" ( dx [ f(p',p")dp'dp"5(x'-x-v p d)5(x"-x-v p ,d).{8l) 
Jx'eA Jx"eB J-oo J 

Now specialize to the case where A is an interval of length I. The total entanglement 
is twice that between A and the real axis to the right of A, which corresponds to 
taking p' < 0, p" > in the above. The integrations over the coordinates then give 
max ((v-pi + Vp")t,£), so that 

/0 />oc 
dp' / dp"f(p',p")(v- p , + Vp„)H(£ - (v- p , + v)*) + 
-oo JO 

/0 /•oo 
dp' / dp"f(p',p") H{{v- P , + v p „)t - I) , (82) 
-oo JO 

where H{x) = 1 if x > and zero otherwise. Now since \v p \ < 1, the second term 
cannot contribute if t < t* = £/2, so that 5^(i) is strictly proportional to t. On the 
other hand as t — > oo, the first term is negligible (this assumes that v p does not vanish 
except at isolated points), and Sa is asymptotically proportional to £, as found earlier. 

However, unless \v\ — 1 everywhere (as is the case for the CFT calculation), 
Sa is not strictly proportional to £ for t > t* . The rate of approach depends on 
the behaviour of f(p' 7 p") in the regions where u_ p / + v p u — > 0. This generally 
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happens at the zone boundary, and, for a non-critical quench, also at p' = 
p" = 0. The exact form of f(p',p") has been exactly calculated only for the 
XY model in a transverse field [124]. The linear increasing followed by (almost) 
saturation has been checked in several lattice models both analytically and numerically 
[113, 53, 125, 126, 127, 128, 129, 124, 130, 131, 132], but we do not have room here 
to discuss the several interesting features that emerged like power-law behaviour for 
large time, periodic time-oscillations etc. 

It is worth mentioning that Eisler and Peschel [130] built a lattice model with 
an exactly linear dispersion relation, and the resulting time-dependent entanglement 
entropy is exactly the one calculated within CFT. It has been also argued that in 
random spin-chains the initial linear growth of the entanglement entropy is replaced 
by a logarithmic one [133]. This is a consequence of the strong scattering among 
quasi-particles and seems to agree with numerical simulations [53]. 

7.1.3. General result for an arbitrary number of intervals. A general result can be 
also derived in the case when A consists of the union of the N intervals (v,2j-\,U2j) 
where 1 < j < N and Uk < Uk+i- Tr//| is given as usual by the ratio Z n /Z n which 
has the form of a correlation function 



in a strip of width 2r . We only need the asymptotic behaviour of this correlation 
function for time t and separations \itj — Uk\ much larger than to- Consequently the 
complicated function T n ,N in Eq. (53) can be set to unity as before. After long algebra 
one arrives to [113] 



If N is finite (or more generally the uj~ are bounded) the second term vanishes for 
sufficiently large t. At shorter times, SA(t) exhibits piecewise linear behaviour in t 
with cusps whenever 2t = Uk — uj, at which the slope changes by ±7tc/6to according 
to whether k — I is even or odd. In the case of an infinite number of regular intervals, 
with Uk — k£, k £ Z, SU(f) exhibits a sawtooth behaviour. 

This behaviour can be explained in terms of the quasi-particles arguments of the 
previous subsection, in which particles entering in and exiting from A entangle and 
disentangle respectively, giving rise to the non-monotonic behaviour. 

The same reasoning applies to other situations, as for example the time 
evolution in the presence of boundaries (particularly relevant for tDMRG that are 
always performed with free boundary conditions). In the simplest instance of the 
entanglement entropy of the segment [0, £] with the rest of the system, the quasi- 
particle argument is easily understood for a perfect reflecting wall at x = 0, for which 
the resulting effective length of the block is 2i and the saturation time t* = £, the 
double of periodic case. This is also easily worked out from the conformal mapping 
z = sin(7rw/2T ) [53, 120]. 

7.2. Local quench 

Suppose we physically cut a spin chain at the boundaries between two subsystems A 
and B, and prepare a state where the individual pieces arc in their respective ground 



(Jj7^(it 2 j-i + in) Y[f n (u 2 j + in)) , 



(83) 



3 j 




(84) 
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Figure 9. Space-time 
region for the density 
matrix for the local 
quench (left) mapped to 
the half-plane (right) by 
means of Eq. (85). z\ = 
ir and 22 = ir + £ 
during the computation 
and in the end r — > it. 
Reprinted with permis- 
sion from [114]. 



states. In this state the two subsystems are unentangled, and its energy differs from 
that of the ground state by only a finite amount. Let us join up the pieces at time 
—t and watch the system evolve up to t = 0. The procedure for the global quenches 
does not apply because the initial state is not translational invariant and will not flow 
under the renormalization group toward a conformally invariant boundary state. 

We can represent the corresponding density matrix in terms of path integral on 
a modified world-sheet. The physical cut corresponds to having a slit parallel to the 
(imaginary) time axis, starting from — oo up to t± = — e — it (the time when the two 
pieces have been joined), and analogously the other term of the density matrix gives 
a slit from T2 = e — it to +oo, like in Eq. (75). Again we introduced the regularization 
factor e. For computational simplicity we will consider the translated geometry with 
two cuts starting at ±ie and operator inserted at imaginary time r. This should be 
considered real during the course of all the computation, and only at the end can be 
analytically continued to it. This plane with the two slits is pictorially represented on 
the left of Fig. 9 where ir corresponds to z\. As shown in the same figure, the z-plane 
is mapped into the half-plane Re w > by means of the conformal mapping 

z I / z\ 2 . w 2 — 1 

w = — h \ / — +1 with inverse z = e . (85) 

e V Ve/ 2w 

On the two slits in the z plane (and so on the imaginary axis in the w one) conformal 

boundary conditions compatible with the initial state must be imposed. 

We consider the time evolution of the entanglement entropy after the local quench 

of two half-chains joined together at the point rry — 0. We consider the four different 

spatial partitions of the system depicted in Fig. 10 among which we calculate the 

entanglement. 



7.2.1. Case I: Entanglement of the two halves. We start with the more natural 
division, considering the entanglement entropy between the two parts in which the 
system was divided before the quench. This is the case when B is the positive real 
axis and A is the negative real axis. Ti p% transforms like a one-point function that 
in the w plane is [2Reu>i]~ d ' 1 . Thus in the z plane at the point z\ = (0,ir) we have 

\ d„ 



{Tn(zi)) 



dw 



dz 



[2Rewi] 



(86) 



that using ewi 



+ \ft 2 — t 2 becomes 
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III) * •--»- ^ 
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Figure 10. The four different bipartitions of the line we consider in [114]. 
Reprinted with permission from [114]. 



Continuing this result to real time r — > it we obtain 

W))=2 n f-^V ■ (88) 



* 2 , 

Using finally the replica trick to find the entanglement entropy we have 



C \og t ^- + ~d 1 . (89) 
6 ae/2 



For f>6we have 



S A (t^>e) = ^\og-+k , (90) 
6 a 

i.e. the leading long time behaviour is only determined by the central charge of the 
theory in analogy with the ground state value for a slit. This could result in a quite 
powerful tool to extract the central charge in time-dependent numerical simulations. 
The constant k is given by k = c[ + (c/6) log(2a/e). 

The behaviour for short time allows instead to fix the regulator e in terms of the 
non- universal constant c^. In fact imposing Sa(£ = 0) = we have e = a/2e~ 6Cl / c . 
Consequently Eq. (89) has no free dynamical parameter, in contrast to the case of the 
global quench. 

This short-time \ogt behaviour has been observed in different situations with a 
local defect [134, 135, 136, 137, 138, 139, 140]. 

7.2.2. Case II: De-centered defect. Let us now consider the entanglement of the 
region r > £ with the rest of the system. In this case Tr p\ is equivalent to the one- 
point function in the plane z at the point z-i = £ + ir as in Fig. 9. Using the conformal 
mapping (9), analytically continuing and taking t,£^> e one finally get [120] 

C, 2£ . 

-log— +c' x t<£, 

( (9i) 

- log — + k t>£, 

b a 

with ko the same as in Eq. (90). The interpretation of this result is direct. Indeed 
at t = the joining procedure produces a quasi-particle excitation at r = 
that propagates freely with the corresponding speed of sound v s that in the CFT 
normalization is v s = 1. This excitation takes a time t — £ to arrive at the border 
between A and B and only at that time will start modifying their entanglement. The 
following evolution for t ^> £ is the same as in Eq. (90). 
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Also the constant value for t < I deserves a comment: it is exactly the value 
known from CFT for the slit in the half- line Eq. (36). This is a non-trivial consistency 
check. Note that a finite e smooths the crossover between the two regimes and makes 
the entanglement entropy a continuous function of the time. 



7.2.3. Cases III and IV: A finite slit. Let us consider again the same physical 
situation as before, but we now calculate the entanglement entropy of A = [0,£] 
and B the remainder. For t < the real negative axis is decoupled from the rest 
and does not contribute to the entanglement entropy that is just the one of a slit in 
half-chain. The entanglement entropy is obtained from the replica trick considering 
the scaling of a two-point function between the endpoints of the slit, that must be 
mapped and analytically continued. After long algebra one gets [120] 

„. ( «(«34)**«0 - „ 

£(M t>£. 



e 2 

Note that compared with Ref. [114] we corrected the behaviour for t < I with the 
generally unknown function T n {x), the boundary analogue of Eq. (42), and we used 
x = 2t/(£ + t) [114] for the four-point ratio. For t > £, only T n (l) — 1 enters and 
the prediction of Ref. [114] remains correct. Using finally the replica trick, we get the 
entanglement entropy {^[{x) = — d n Fn{x)\ n =i) 



S A = < 



c , t c, £ c , £-t ~, f 2t \ „ . 

-ln- + -ln- + - In 4 h T\ + 2c', t < 

3 a 6 e 6 (, + t 1 \£ + 1 J 1 

c I 

- In - + 2c[ t > 

3 a 



(93) 



The crossover time t* = I is again in agreement with the quasi-particles interpretation. 

There are several interesting features of this result. For very short time t <SL £'\t 
reduces to the £ — oo case Eq. (90) as it should. The leading term for t > I is just the 
ground state value for a slit in an infinite line. However the subleading term is not 
the same, signaling that for long time the system still remembers something of the 
initial configuration as a boundary term that is unable to "dissipate" . Since the extra 
energy never dissipates under unitary evolution, there is no reason for the constant 
terms to be the same. According to Eq. (38) these two constant terms are the same 
only when .9 = 1, as it is the case for the Ising model with free boundary conditions. 

Another interesting feature is the behaviour for t < £. This is very similar to the 
form proposed in Ref. [135] to fit the numerical data, i.e. 

S A = C flog£+ C j \og(t/£) + | log(l - t/£) + k' . (94) 

Only the terms in t + I and T' x were missing in Ref. [135]. However these behaves 
smoothly for < t < I and their effect can be well approximate by a constant factor 
in k'. The results of the fit are Co — 1 + C2, c\ ~ 1, and c 2 ~ 1/2 that are exactly our 
predictions for c = 1 . 

For the most general case of a slit A — \l<2,l{\, we remand to the original paper 
[114]. We stress here that the result for short times (Eqs. (35) and (36) there) are 
generically incorrect because we assumed F n (x) — 1, that we learn successively not 
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to be the case. The result for t > £\ are instead correct, because in the regime only 
T(V) = 1 enters and in this case we have 

S A (t>e 1 ) = ^ln e -^—^+2c' 1 , (95) 

o (X 

that is the same as for case III. The correct results for t < £\ in the various regimes, 
can be read from Eqs. (48), (50) and (51) in Ref. [114] from the general results of 
the two-point functions of generic primary operators. We mention that we could use 
non-equilibrium calculations/simulations to determine the unknown function T n {x). 
It is still not clear if this could be effective. 

An interesting generalization of the local quantum quench in models with 
gradients has been provided in Ref. [141], where again all the entanglement evolution 
can be understood in terms of the quasi-particle picture. 



7.2.4- Decoupled finite interval. A natural question arising is how the results we 
just derived change when we introduce more than one defect in the line. It is 
straightforward to have a path integral for the density matrix: we only need to have 
pairs of slits for — oo to — ie and from ie to +ioo everywhere there is a defect. However 
it becomes prohibitively difficult to treat this case analytically. In order to begin to 
understand the case when a finite interval interval is initially decoupled, we consider 
the case when it lies at the end of a half-line. 

So, let us consider a semi-infinite chain in which the A subsystem is the finite 
segment (— £, 0) and the B is the complement (0, oo) and with the initial defect 
at rrj = 0. The space-time geometry describing this situation is like the one just 
considered, with a wall at — £ + iy (y real) that represents the boundary condition. 

In these circumstances the inverse conformal mapping between the z plane and 
the half-plane can be worked out using the Schwarz-Christoffel formula. After long 
algebra one obtains [114] 



z(w) = i 



■ log(iw) + b- 



1 



1 



(96) 



with the parameter b related to £ and e in a non-algebraic way. (A slit in the full line is 
closely related to this transformation, the last piece is replaced by (w 2 — l)/(w 2 + 1).) 
Unfortunately the mapping (96) is not analytically invertible and its exact use is 
limited to numerical calculations that do not help us, since we need to perform an 
analytical continuation. However, even if not completely justified, we can take the 
limit £ e, before the analytical continuation to real time obtaining [114] 

1 



w = — i exp 



me 
~2£ 



2 /e 2 + 1 - z/e) 



.£, ,. \ .e 2 n 
z = i— log(iw) + l- 

7T 



I 4 log(iw) ' ^ 

This approximation is expected to fail for t > £. It is easy to perform the mapping, 
continuing to real time r — > it, and for t >• e we have [114] 



(T n (ir)) = c r , 



(98) 



_4£ isin(7rf/2f)_ 

Clearly this cannot make sense when the argument of the power-law becomes negative 
(i.e. for t > 2£), signaling the expected failure of Eq. (97). The time when this 
approximation fails cannot be understood from this calculation, but only in the 
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comparison with explicit results in real time (or by the exact use of Eq. (96)). Using 
the replica trick, for the entanglement entropy we obtain 

S^-logl-tsin-j-fSi. (99) 

One is tempted to assume that this result can be correct for t < I and that for 
larger time it saturates as suggested by the quasiparticle interpretation. In the case 
of an initially decoupled slit of length I in an infinite chain Eq. (99) is still valid with 
the replacements 21 — > £ and c/6 — > c/3 as follows from a simple analysis [136]. The 
validity of this equation has been carefully tested for the XX chain in Ref. [136] , finding 
very good agreement for all t < £, confirming the naive expectation. In Ref. [47] a 
more complicated kind of defects has been investigated, and the results always agree 
with Eq. (99) when describing a conformal hamiltonian. 

8. Local quench, quantum noise and measuring the entanglement 

The entanglement entropy has been revealed to be a useful quantity for a deep 
theoretical understanding of extended quantum systems, especially in connection 
with criticality and topological order (see the review by Fradkin in this volume [41]). 
However, the final check and goal of any theory is the comparison with experiments. 
The intrinsic non-local nature of the entanglement entropy makes any attempt at an 
experimental measurement difficult, if not impossible. Some bounds relating Sa to 
thermodynamic observables have been derived [142], but this is still far from being an 
operational measure. 

However, there is one recent interesting proposal to measure the entanglement 
entropy out of equilibrium, in the setup of the local quench we have just described 
[7, 143, 144]. The main idea of Klich and Levitov is to relate the entanglement between 
two half-chains to the distribution of the electrons passing towards the contact between 
them. They considered (as we did above) two semi-infinite chains (which are leads in 
actual experiments) initially disconnected and then at some time to joined together, 
allowing the passage of electrons (if the leads are two Fermi seas, the quasiparticles of 
above are real electrons). The transport at this quantum point contact is described by 
the theory of quantum noise. This approach describes the probability distribution of 
transmitted charge using the generating function x(A) = X^-oo Pne tXn , where P n 
is the probability to transmit n charges in total. The function x(A) can be written in 
terms of cumulants C m 

l 0gx(A) = g (100) 

( m! 

m— 1 

The fundamental point is that the constants C n are measurable quantities (C 2 is 
measured in routine experiments, and also C n up to n = 5 have been measured in 
more difficult experiments). 

The main result of Ref. [7] is to establish a relation between the cumulants C n 
and the entanglement entropy of the two-halves 

c V a ™r n f ( 2w ) m \ B ml meven 

S * ~ 2. a ™ = { 0, m odd ' (101) 

m>0 

where B m are Bernoulli numbers. For quantum noise generated in the contact 
switching on (at to) and off (at ti), the current fluctuations are gaussian (C m7 t2 = 0), 
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with a variance C 2 = -\ log tl ~ t ° , where r is a short time cutoff set by the contact 
switching rapidity. Combined with Eq. (101) this gives entropy Sa ~ (1/3) log \ t-\_— to I ■ 
In Ref. [7] this has been put in direct relation with the standard formula Sa = c/3 log^ 
(c = 1 of free electrons). However, we have seen in the previous section Eq. (90), that 
c/3 log t is a key feature of the local quench that comes from the specific time-scale 
e whose analogous here is r. In Ref. [7] also the reaction of the system to a train of 
pulses (i.e. periodic switching on and off of the contact) has been considered. When 
two point-contacts are activated at the same time the response should be given by 
Eq. (99). 

We finally stress that the previous analysis is valid for free electrons, and it is 
unclear at present how the treatment must be properly modified in general to take into 
account interactions to describe other universality classes. A first calculation for the 
Luttinger liquid theory of the quantum Hall point contact showed that the measured 
noise is always logarithmic, with a prefactor not given by the central charge, but by a 
filling v dependent constant [140] 



X(A) = exp 



A 2 v At 



(102) 



for At = ti — t 3> t. Furthermore, it has been shown that for the Ising model 
the noise at a point contact is algebraic instead of logarithmic [140], suggesting that 
the relation between the full counting statistics and the entanglement could not be 
structural. 
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